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Abstract 

We present two simple analytical methods for computing the gravity wave 



^ ' contribution to the cosmic background radiation (CBR) anisotropy in infla- 

tionary models; one method uses a time- dependent transfer function, the 

^ , other method uses an approximate gravity- wave mode function which is a sim- 

H 

^ ' pie combination of the lowest order spherical Bessel functions. We compare 

the CBR anisotropy tensor multipole spectrum computed using our methods 
^ ' with the previous result of the highly accurate numerical method, the "Boltz- 

mann" method. Our time-dependent transfer function is more accurate than 
the time-independent transfer function found by Turner, White, and Lidsey; 
however, we find that the transfer function method is only good for / < 120. 
Using our approximate gravity-wave mode function, we obtain much better 
accuracy; the tensor multipole spectrum we find differs by less than 2% for 
/ < 50, less than 10% for / < 120, and less than 20% for / < 300 from 
the "Boltzmann" result. Our approximate graviton mode function should be 
quite useful in studying tensor perturbations from inflationary models. 



PACS index numbers: 04.30.+X, 98.80.Cq, 98.80.Es 
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1. Introduction 



The cosmic background radiation (CBR) anisotropy places stringent constraints on the- 
ories of the early Universe. Among these theories, the best studied are the inflationary 
models, which are strongly motivated because they solve the famous problems (the flat- 
ness problem, the smoothness problem, the structure formation problem) of the standard 
cosmology. Tensor (gravity-wave) and scalar (density) metric perturbations are generated 
in the very early Universe due to quantum fluctuations arising during inflation. Both 
tensor and scalar perturbations contribute to anisotropy in the temperature of the CBR. 
While the scalar contribution to the CBR anisotropy involves more complicated physics, 
the tensor contribution to the CBR anisotropy arises only through the Sachs- Wolfe effect 
|l|] as follows. As photons of the CBR propagate toward us from the last scattering sur- 
face, their paths are perturbed by the metric perturbations due to the primordial gravity 
waves. The perturbed energies of these photons result in temperature fluctuations in the 
sky that we observe. The CBR temperature fluctuation is conventionally expanded into 
spherical harmonics: 



In this paper, we present two simple analytic methods of computing the tensor contribu- 
tion to the variance in the CBR temperature multipole moments, (|a/mP); one method 
makes use of a time-dependent transfer function, the other uses an approximate gravity- 
wave mode function which is a simple combination of the lowest order spherical Bessel 
functions. Our methods provide adequate accuracy for normalizing the tensor perturba- 
tions arising from inflationary theories to the observable CBR anisotropy. Our method 
using the approximate gravity-wave mode function is accurate enough to be used in 
studying the CBR anisotropy tensor multipole spectrum to large / (/ < 300). 
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(1.1) 
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We use the results by Dodelson, Knox, and Kolb as the standard for compari- 
son. They considered a Universe with both matter and radiation, and used numerical 
methods to evolve the photon distribution function using first-order perturbation the- 
ory of the general relativistic Boltzmann equation for radiative transfer. Their results 
(which we refer to as "Boltzmann") should be equivalent to the results obtained using 
the Sachs- Wolfe formula (see below). The "Boltzmann" method has no simple analytical 
formulation, and the exact result using the Sachs- Wolfe formula involves complicated 
spheroidal wavefunctions 0]; hence it is of great interest to find a simple transfer func- 
tion, or a simple approximate graviton mode function, which can be used analytically to 
compute the CBR anisotropy tensor multipole moments to sufficient accuracy. 

Our first method has been motivated by the tensor transfer function found by Turner, 
White, and Lidsey (TWL) [Q. Their transfer function takes into account the effect of the 
Universe becoming matter-dominated gradually. The tensor multipole moments obtained 
using their transfer function, however, differ from the "Boltzmann" results by up to over 
30% for / < 100. The reason for this substantial discrepancy is that their transfer function 
contains no time evolution, although they did give the expression for the time-dependence 
of the transfer function for short wavelength modes. We modify their transfer function 
by accounting for the difference in time evolution between long and short wavelength 
modes. 

Our second method has been motivated by intuition. Since the gravity-wave mode 
function is analytically known for both matter and radiation dominated eras, it should 
be possible to construct a simple approximate gravity-wave mode function by smooth 
interpolation. As expected, the resultant mode function is much closer to the true mode 
function than the much used matter-dominated mode function. When used in computing 
the CBR anisotropy tensor spectrum, the approximate gravity-wave mode function gives 
very accurate results. 
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2. Gravity wave contribution to CBR anisotropy 

The rms temperature fluctuation averaged over the sky for a given experiment is given 
by 

((f)^)-g^(l.™P)«'. (2.1) 

where Wi is the appropriate response function for the experiment. For an experiment 
with two antennas of Gaussian beam width a separated by angle 6, the temperature 
difference between the two antennas is measured; the response function is 

Wi = 2[1 - P,(cos^)] e-('+i/2)V^ (2.2) 

We have followed the notation of Ref . . 

Tensor perturbations generated by inflation are stochastic in nature Let us ex- 
pand the gravity-wave perturbation in plane waves 

/.,,(x, r) = (27r)-3 J d'k hUT)e^, e-''^'^ (2.3) 

where e*;, is the polarization tensor and i = x, + in the transverse traceless gauge (in 
which /loo = hoj = 0). We have 

hiiO) = A{k)a'ik), (2.4) 
where a*(k) is a random variable with statistical expectation value 

(a,(k) a,(q)) = k'^'S^'^Xk - q) (2.5) 
and the spectrum of gravity waves generated by inflation is 

where V is the value of the inflaton potential when the mode with comoving wavenumber 
k crosses outside the horizon during inflation. 



The tensor contribution to the variance of the multipoles is given by 
(|a,„|2) = 367r2|±|l| J dk A\k) \Fi{k)\\ 

where 

jiik{Tnaw - r)) 



(2.7) 



1 /"Tnow Q 

^^^^^ ^/k Xlss dr 



1 hlir) 
3 hliO) 



(2.8) 



[k{Tnow - r)f ' 

where Tnow is the conformal time today, tlss is the conformal time at last scattering. We 
need to find the graviton mode function h1_{T). 

Infiation gives rise to a spatially fiat and perturbed Friedmann-Robertson- Walker 
universe with the metric 



(2.9) 



where rj^i, = diag{l, —1, —1, —1), h^y is a small perturbation, and r is the conformal 
time. The cosmic scale factor R{t) is 



i-cq 



(2.10) 



for a Universe with both matter and radiation. We have defined tq = 2Hq 1 + -Req 
At matter-radiation equality, -R(req) = i?eq = 4.18 x lO^^/i^^, and Teq/ro = [V2 — l]Rl/'^ 



eq 



Today R^o^ = 1, T^ow/ro = + Rcq - ^/R^- At last scattering, i?Lss = V(l + ^^lss), 
tlss/tq = 



Ro 



Req- 



l-LSS I -'•■eq 

The gravity-wave perturbation satisfies the massless Klein-Gordon equation 



K + 2 



hi + h 



2 ui 
k 



0, 



(2.11) 



where fc^ = k ■ k and the overdots denote derivatives with respect to r. 

A gravity- wave mode with wavenumber k crosses inside the horizon at /cr ~ 1. Before 
it crosses inside the horizon, fcr <^ 1. Eq. (p.ll| ) gives us hl^{T) = for /cr ^ 1, i.e., the 
gravity-wave mode is frozen before horizon-crossing. We can take /^^(O) = as the 



initial condition for Eg. ( |2.11| ). For modes that cross inside the horizon during radiation 



dominated era (r <^ Teq, R{r) = 2^J Rcc^t /tq), the exact solution is /^^(t) = ^k(O) io(^''"); 
for modes that cross inside the horizon during matter dominated era (r ^ Tgq, R{t) = 
{t/tqY) the exact solution is hl-ir) = /^^^.(O) 3ji(A;r)/(A;r). Here jo(^) = sinz/z and 
ji{z) = sinz/z"^ — cos z/ z are spherical Bessel functions of order zero and one respectively. 

3. First method: time-dependent transfer function 

Since the contributions to the tensor multipole moments are dominated by gravity waves 
which have entered the horizon recently 0, let us write 0] 

"3ji(A;r) 



(3.1) 



kr 

where Tr{k/keq) is the amplitude transfer function which accounts for the effect of short- 

I'eq — 'cq 



wavelength modes entering the horizon during radiation-dominated era, and ken = r^,, ^ 



Eq.(p^.8|) becomes 

F.(k) = -k^' r-drri;W ^*g^Ww-;)' . (3.2) 

Note that there is no dTT-{k)/dT term in the above expression, because dh\^{T)/dT and 
h\Jj^) are related to 3ji(A;r) / [kr) and 9[3ji(A;r) / (A;r)]/9r by the same amplitude transfer 
function Trik). 

The transfer function at time r can be found by integrating Eg. ( |2.1lD numerically 
from r = to r. Today's transfer function is [^] 

T.o(l/) = To(i/) = [1.0 + 1.3% + 2.5y2]'/', (3.3) 
where y = k/kcq. 

Since the Universe became matter-dominated gradually, the transfer function in 
Eg. ( ^.11 ) should obviously depend on time. Once a mode is well inside the horizon 
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{kr ^ 1), /^^(t) oc cos{k:r)lR (see Eq. Since 3ji{kT)/{kr) is the exact mode 
function for R{t) = (r/ro)^, the transfer function for modes with k ^ k^q at an early 
time r is given by [Q 



i?(r) 



eq- 



(3.4) 



The above formula is in very good agreement with numerical results for k > kgq and 
T-Lss <r <To. 

On the other hand, modes with k -C fceq entered the horizon during matter-dominated 
era; the transfer function for these modes should have negligible time dependence. Let 
us write 



Tr{k/kcq) = To{k/kcq) Ti(r, k/k 



(3.5) 



where To^k/keq) is given by Eq. (p^) , and Ti{t, k/k^q) can be written as 
Ti(r, fc/fceq) = w{k/k,q)A{T) + [1 - W{k/k,q)]. 



(3.6) 



where A{t) = {t/tqY /R{t). w{k/kcq) for /c <^ /cgq, and w{k/kcq) ^ 1 for A; fc, 



eq- 



The simplest choice is 



t(7(/i;//i;eq) = 1 — exp —^{k/k, 



(3.7) 



where 7 and A are constants. 

Since we use the graviton mode function from matter-dominated era, it is consistent 



to use tlss — tq-vZ-Rlss (which is the matter-dominated limit of the correct expression) 
as the lower limit of integration in Eq. (|2.8|) . With 7 = 0.9 and A = 0.45 in Eq. (|3.7| ), the 
multipole moments computed using our transfer function agrees reasonably well with the 
"Boltzmann" result for / < 50 (see Fig. 1). To get better result at larger /, we can use a 
slightly complicated weight function 



w{k/keq) = (l-exp -7(/c//Ccq)^ ) (l - cxp -ndk/k^q 



(3.8) 



where ric and Uc are constants. With 7 = 0.9, A = 0.45, ric = 20, i/c = 0.7, the 
multipole moments computed using our transfer function agrees reasonably well with the 
"Boltzmann" result for / < 100. In Figure 1, we plot the CBR anisotropy tensor multipole 
spectrum computed using the TWL transfer function (dotted line), the transfer functions 
with weight functions given by Eq.( p.7|) (dashed line), and Eq.( p.8|) (dot-dashed line). The 
solid line is the "Boltzmann" result. 

It can be argued that using the matter-dominated limit for tlss distorts the ionization 

is 



history of the Universe, since the correct expression tlss = tq \/ Rlss + -Req — y-Rcq 



smaller than the matter-dominated limit tlss — tqV-Rlss by a factor of 2/3. However, 
using the correct expression for tlss increases the multipole moments for an amount which 
increases from 10% at / = 20 to 300% at / = 100. The reason for this dramatic effect is 
that the graviton mode function we use is an extremely bad approximation to the true 
mode function at the era of last scattering. The contribution to the multipole moment 
for a given / is dominated by the wavenumber at which Fi{k) (see Eq.( |2.(j| )) peaks; since 
Fi{k) peaks at kro ~ /, i.e, k/keq ~ 2.678/i^^ / x 10^'^ larger / multipole moments are 
dominated by the contribution from larger-wavenumber graviton modes, which entered 
the horizon at earlier times, when the true graviton mode function deviates greatly from 
the matter-dominated graviton mode function that we use. Our transfer function can 
not correct for this effect even with time dependence included, because at early times 
(around the last scattering) the graviton mode function has only a smaller number of 
oscillations in k, while our transfer function only accounts for the difference in average 
amplitude between the true and the matter-dominated mode functions. 

By using the matter-dominated limit for tlss? we are effectively truncating the inte- 
gral over conformal time r in the expression for the multipole moment; it is not surprising 
that this enables us to get multipole moments (for / < 120) which are not far off from 
the "Boltzmann" results, since we are cutting off the r integral at small r where the 
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matter-dominated mode function is most inaccurate. The multipole moments we obtain 
by using the matter-dominated hmit for tlss are therefore physically consistent. 

The transfer function method is limited to / < 120. For larger /, the phase difference 
between the matter-dominated graviton mode function and the true mode function be- 
comes important, which results in significant deviation between the multipole spectrum 
computed using any transfer function and the multipole spectrum from the "Boltzmann" 
method. 

4. Second method: approximate gravity- wave mode 
function 



If we want to use the correct expression for tlss, we must use a new mode function 
which better approximates the true graviton mode function at small r than the matter 
dominated mode function 3ji(A;r)/(A;r). 

Obviously, one can construct an approximate solution which interpolates smoothly 
between jo^kr) for r ^ r^q and 3ji(A;r)/(A;r) for r ^ Teq. Let us write 



d_ 



[l-w{T)]T^{k/k, 



d 



0, 



-ki[i-w{T)mk/k,q) 



kr 
kr 



+ voir 



djoikr) 
dr ' 

+ w{T)ji{kT) 



(4.1) 



where w(r) ^ for r > r^q, and w{t) ^ 1 for r < Teq. T^k/k^q) = [{t / tq^ / R{t)] To{k/k 
(see Eq.(|3.4|)), it correctly accounts for the difference in average amplitude between the 
matter-dominated mode function and the true mode function. 

Since joikr) accurately gives both the amplitude and the phase of the true mode 
function at large k, we should "turn off" Tr^Xk/k^q) (which accounts for the amplitude 
difference between 3ji(A;r)/(A;r) and the true mode function) at large k. We can replace 
Eq.(O) with 



eq/ 



-ay" 



1.0 + 1.34y + 2.5y2]^^^ + 



-ay- 



y = k/k 



eq- 



(4.2) 



To compute the multipole moments, we make the following substitution in Eq. (|3.2|) : 

^r{k) ^ [1 - w{t)] (fc/fccq) — ^ + W[t) — (4.3) 

A good choice for w{t) is 

ti;(r) = e-"("/"^'')^ (4.4) 

where a and (3 are constants. 

Fitting the resultant tensor multipole spectrum to that found by the "Boltzmann" 
method, we find a = 0.2, (3 = 0.65, a = b = A (the fitting is more sensitive to a and 
f3 than to a and b). The agreement between our result with the "Boltzmann" result is 
impressive: better than 2% for / < 50, better than 10% for / < 120, and better than 
20% for / < 300. In Figure 2, we plot the CBR anisotropy tensor multipole spectrum 
computed using the matter-dominated graviton mode function with the TWL transfer 
function (dotted line); our approximate graviton mode function (given by Eqs.( [4.1[ ), ( |4.2| ), 
and ( |4.4| )) (dashed line). The solid line is the "Boltzmann" result. 

The multipole spectrum computed using the "Boltzmann" method has a peak at 
/ ~ 217 (see Figure 2). The multipole spectrum computed using our approximate mode 
function has a peak at / ~ 213, while the multipole spectrum computed using the matter- 
dominated mode function (with or without transfer function) has a peak at / ^ 180 (see 
Figure 2). Our approximate mode function gives much more accurate phase information 
than the matter-dominated mode function. 

The successful application of our approximate graviton mode function in computing 
the tensor multipole spectrum stems from the fact that it is rather close to the true 
graviton mode function. In Figures 3 and 4, we plot the conformal time derivatives of 
the true mode function (solid line), our approximate mode function (given by Eqs.( [4.1| ), 
( [4.2| ), and ( [4.4D ) (dashed line), and the matter dominated mode function (dotted line). 
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Figure 3 shows that at r = tlss? our approximate mode function is much closer to the 
true mode function than the matter-dominated mode function, for all wavenumbers k. 
Figure 4 shows that for k = k^q = t~^, our approximate mode function is also much 
closer to the true mode function than the matter-dominated mode function. 

5. Generalization 

The CBR anisotropy tensor multipole spectrums in Figures 1 and 2 are for the standard 
values of h = 0.5, = 0.05, Qq = 1; and for a scale-invariant primordial spectrum of 
gravity waves. It is straightforward to use our methods to compute the tensor multi- 
pole spectrum for a non-scale-invariant primordial spectrum of gravity waves, just use 
the corresponding A{k) (see Eq.( p.6|) , which is constant in the scale-invariant case) in 
Eq.Q. 

Next, we consider the cases when some or all of the parameters h, Q-q, and Qq are 
different from the standard values {h = 0.5, Qb = 0.05, Qq = 1). 2;lss is given by 

l + ZLss^llOO(fio/^^Br''. (5.1) 

As long as f^o = 1; the cosmic scale factor is given by Eq. (|2.10| ); our previous formalism 
applies with zi^ss given above. In Figure 5, we plot the CBR anisotropy tensor multipole 
spectrum computed using our approximate graviton mode function (given by Eqs. ([4.1|) , 



(13), and (0)) for = 1, ^^b = 0.05, h = 0.8 (sohd line) and h = 0.5 (dashed line). 
If f2 = fio + = 1, but Qq < 1, we have 
r 1 dR 



^0 2 Jo jR^^ + R+(njno)il + Req) i?^ 



i?eq = 4.18 X 10"^(fio/i^^~^ 

L — -1 _ "^0 ''^ 

"^eq — ''"cq 



(V^-I)^i?eq 
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Note that the conformal time today Xnow 7^ tq- For h = 0.8, Qq = 0.2, Q\ = 0.8, Tjiow/to — 
0.85. A sizable cosmological constant term significantly alters the late time evolution of 
the cosmic scale factor, hence 3ji(A;r)/(A;r) is no longer the mode function at late times. 
However, our approximate graviton mode function is still a good approximation of the 
true mode function in the context of computing tensor multipole spectra, which are not 
very sensitive to the late time graviton mode function |Q. In Figure 6, we plot the CBR 
anisotropy tensor multipole spectrum computed using our approximate graviton mode 
function for Qq = 0.2, = 0.8, = 0.02, h = 0.8 (solid line), and Qq = 1, ^^b = 0.05, 
h = 0.5 (dashed line). 

We expect the tensor multipole spectrum computed using our approximate graviton 
mode function (given by Eqs.(|3), G), and (U), with a = 0.2, /3 = 0.65, a = b = A) 
to have fairly good accuracy when the cosmological parameters h, Qb, and Qq differ from 
the standard values of h = 0.5, Qb = 0.05, Qq = 1. 

6. Discussion 

We have presented two simple and straightforward analytical methods for computing the 
CBR anisotropy tensor multipole spectrum; one method uses a time-dependent transfer 
function, the other method uses an approximate gravity-wave mode function which is 
a simple combination of the lowest order spherical Bessel functions (given by Eqs. ([4.1|) , 
(13), and (|]3), with a = 0.2, (3 = 0.65, a = b = A). Both methods give much better 
accuracy than using the matter-dominated mode function with the time independent 
transfer function of Turner, White, and Lidsey [Q. Our approximate graviton mode 
function method is especially promising, since it gives a tensor multipole spectrum which 
is impressively accurate (difference with the "Boltzmann result" is less than 2% for / < 50, 



12 



less than 10% for / < 120, and less than 20% for I < 300). 

Many attempts have been made in the past in finding a good approximate graviton 
mode function, for instance, via the sudden approximation (assuming the transition 
from radiation-domination to matter-domination to be instantaneous) In addition, 
Ng and Speliotopoulos explored the possibility of finding a good approximate graviton 
mode function via the WKB approximation [||; however, they did not use the mode 
function they found in computing the multipole moments, they used the mode function 
found by numerical integration instead. Furthermore, Koranda and Allen found the 
exact graviton mode function in terms of spheroidal wavefunctions Q]. The advantage 
of our approximate graviton mode function is that it takes into account the gradual 
transition from radiation-domination to matter-domination in the Universe; it can be 
used to compute the CBR anisotropy tensor multipole spectrum rather accurately; and 
it involves only the first and second order spherical Bessel functions. 

In a different direction of the CBR anisotropy tensor multipole spectrum calculation, 
the "Boltzmann method" evolves the photon distribution function numerically, using 
first-order perturbation theory of the general relativistic Boltzmann equation for radiative 
transfer 0. We have used the CBR anisotropy tensor multipole spectrum from the 
"Boltzmann method" as the standard for comparison, but we do not think that it is 
practical to use the numerically rather involved "Boltzmann method" in all cases. When 
we study the tensor perturbations from a large number of inflationary models, it is much 
more convenient to use the Sachs- Wolfe formula with a simple yet reasonably accurate 
graviton mode function, such as the one presented in this paper. 
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Figure Captions 



Figure 1. The CBR anisotropy tensor multipole spectrum computed using various 
transfer functions, compared to the Boltzmann result (solid line). 

Figure 2. The CBR anisotropy tensor multipole spectrum computed using our ap- 
proximate graviton mode function (given by Eqs.( [4.1|) , ( [4.2|) , and ( [4.4| )), compared to the 
Boltzmann result (solid line). 

Figure 3. The conformal-time derivatives at r = tlssj of the true graviton mode 



function (solid line), our approximate mode function (given by Eqs.([4.1|), (^.2|), and 
( [4.4D ) (dashed line), and the matter-dominated mode function (dotted line). 

Figure 4. The conformal-time derivatives with k = keq = t^, of the true graviton 
mode function (solid line), our approximate mode function (given by Eqs. (|4.1| ) 



and (|4.4| )) (dashed line), and the matter-dominated mode function (dotted line). 

Figure 5. The CBR anisotropy tensor multipole spectrum computed using our ap- 
proximate graviton mode function (given by Eqs. (|4.1| ), ( [4.2|) , and (|4.4| )) for f2o = 1, 

= 0.05, h = 0.8 (sohd line) and h = 0.5 (dashed line). 

Figure 6. The CBR anisotropy tensor multipole spectrum computed using our ap- 
proximate graviton mode function (given by Eqs.( [4.1| ), ( |4.2| ), and ( [4.4[ )) for = 0.2, 
Qa = 0.8, Qb = 0.02, h = 0.8 (solid line) and Qq = 1, fig = 0.05, h = 0.5 (dashed line). 
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